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^AJ , This is a short presentation of our talks given at eRHIC Workshop at the BNL. We give 

Qs ' here a status report of our attempts to understand how gluon saturation will manifest itself 

in deep inelastic scattering with nuclei. This summary reflects our current understanding 

and shows directions of our research rather then a final answer to the question. Nevertheless, 

we are able to share with our reader our tentative answer to the question: "Why do we need 

^ , to measure DIS with nuclei and why these data will be complementary to the information 

r^ ' obtained from proton DIS". 
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I. INTRODUCTION: WHAT ARE THE SCALES IN PHOTON-NUCLEUS DIS? 



(N 

o 

^D I The main goal of these notes is to examine if and how we can observe the phenomenon of gluon saturation 

r~| ' in DIS with nucleus. We present here only a status report of our attempts to clarify this subject, which is far 

^ i' away from being complete. It, rather, indicates the directions of our searches. Much more work is needed to 

develop a reliable approach so as to finalize our recommendations concerning experiments the most sensitive 

to the gluon saturation. 
f-* ', We start with the general approach to photon-nucleus interaction, developed by Gribov [y who suggested 

following time sequence of this process: 



1. First, the 7* fluctuates into a hadron ( quark - antiquark ) system well before the interaction with 
the target; 

2. Then the converted quark- antiquark pair ( or hadron system) interacts with the target. 
Generally, these two stages result in the following formula for the cross section 

atot{-f*+A) = Y^ I'i'nl^ atot{n + A]x) , (1) 

n 

where ^„ is the wave function of the system, produced in the first stage of the process. 

A. Separation scale r^'' ~ 1/Mo 

This scale is a typical distance which separates the pQCD approach from the non-perturbative one. 
Roughly speaking, for shorter distances than r''^^ , the QCD running coupling constant can be considered 
as a small parameter while for longer distances cts{r±) is large and we cannot use the powerful methods of 
pQCD. Table 1 demonstrates how this scale works in our particular model to incorporate the long distance 
physics 



Table 1 



Perturbative QCD 
short distances 

r±_ < 

DOF: colour dipoles H 
'^„: QED for virtual photon 

atot{n,x) = a{r1,x) 

Glauber- Mueller Eikonal for (j{r\,x) 



Sep 



non-perturbative QCD 
long distances 

< r±_ 

DOF: constituent quarks H 

*&„: generalized VDM for hadronic system 

atot{n, x) = a{qq -^ qq; x) 

Regge phenomenology for a^q + q ^> q -\- q; x) 



It is important to notice that the separation scale mostly relates to the produced hadronic {qq) ) system 
and does not depend on the properties of the target ( in particular, the atomic number ). From Table 1 one 
can write for short distances ( r^ < r 



Sep 



<ytot{l*p) = I d r± I dz|vI/(Q ■,r^,z)\ atot{r^,x). 



(2) 



B. Saturation scale r^ 



l/Q,(x;A) 



At low X and at r± < r'^^^ we believe 



that The system of partons always passes the stage of hdQCD 
( at shorter distances ) before it goes to the black box, which we call non-perturbative QCD, and which, in 
practice, we describe in old fashion Reggeon phenomenology. At the hdQCD stage we have to observe a 
parton system with sufficiently small typical distances ( r^°* w l/Qs[x\A) ) at which the QCD coupling 
constant is still small {as{rj^*) ^ 1), but the density of partons is so large that we cannot use the pQCD 
methods in our calculations. The picture of the parton distribution in the transverse plane is shown in Fig. 
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The estimate of the value for the saturation scale is obtained ||g 

fsi^)) 



from the equation 



K 



i-K^asA xG(x 

X 



2g2(a;) 



nR\ 



(3) 



where A and Ra are the atomic number and radius of the nucleus. Eq. (13) has a simple physical meaning 
giving the probability of the interaction between two partons in the parton cascade. Namely, such an 
interaction will stop the increase of the parton density due to parton emission, which is included in the 
DGLAP evolution equations [||. 

It is important to notice that the saturation scale strongly depends on A {Qs{x; A) ex Ae [|l^ -^ A^ [1^ ). 



C. The theory status 

In cA deep inelastic scattering we want to find the high density parton system which is a non-perturbative 
system but which can be treated theoretically. It should be stressed that the theory of hdQCD is in a very 
good shape now. Two approaches have been developed for hdQCD: the first one |l2| is based on pQCD ( see 
GLR and Mueller and Qiu in Ref. gl) and on the dipole degrees of freedom ||], while the second [Q uses 
the effective Lagrangian, suggested by McLerran and Venugopalan M . As a result of this intensive work we 
know now the nonlinear equation which governs the QCD evolution in the hdQCD region |]1J]. We have not 
developed simple methods to estimate an effect of hdQCD on the experimental observables and have to use 
a model approximation, but we want to emphasize that this is a temporary stage of our theory which will 
be overcome soon. 



II. HERA: RESULTS AND PUZZLES. 

We start answering the question: "why do we need a nuclear target to find the hdQCD phase" with a 
summary of what we have learned from HERA. 
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A 



K=l 




K« 1 



000 



> Q' 



FIG. 1. Tfie parton distribution in the transverse plane. The curve shows the saturation scale Qs{x;A) 

F2 - the most striking and significant result from HERA is the increase of F2 at low x |Q. The 
interpretation of the F2 data in terms of the DGLAP evolution equations leads to sufficient large value 
and a steep behaviour of the gluon structure function at low x. xG{x,Q^) turns out to be so large 
that K, our new order parameter, exceeds unity in a significant part of the accessible phase space ( see 
Fig. I). 




FIG. 2. 



Diffractive production - three important results have been observed at HERA: (i) the diffractive 
production gives a substantial part of the total cross section, about 10 -^ 15 % at Q^ « lOGeV^; (ii) 
the energy behaviour of the diffractive cross section has an intercept larger that the intercept of the 
soft Pomeron, namely, (Tdiff oc (1/a;)^'^^ with Ap > Asoftp ~ 0.1 [|6|, and (iii) the ratio adiff /(Jtot is 
a constant versus energy in HERA kinematic region. From Fig. |l|, one can see that a hadron looks as 
a diffractive grid with typical size of the order of r""*. Therefore, we expect that diffractive processes 
originate from a rather small distances. This fact leads to a natural explanation of the energy behaviour 
of the diffractive cross section. 

Matching between soft and hard processes. The experimental data on 7*p cross section at 
small Q^ allows to test different models for the matching of the soft and hard interactions. 

The dedicated beautiful measurement of the F2 slope {dF2/ hiQ^) gives us a hope to find the 
saturation scale by observing the movement of the maxima in Q^ - behaviour at fixed x. The exper- 
imental data show a considerable deviation fron the DGLAP analysis at Q^ < 1 ^ 3GeV^. However, 
the current data can be described in two different ways, either due to a gluon saturation or due to a 
probable matching between soft and hard at rather large momenta ( about 1-2 GeV ) H. It should 
be noticed, however, that the J/^ production can be easily described taking into account shadowing 
corrections confirming a gluon saturation hypothesis M. 



We listed above the most important HERA observations which indicate a possible saturation effect. To 
illustrate this fact we will demonstrate that a simple parameterization of Golec-Biernat and Wusthoff fT^, 
which includes the saturation, works well. They found an elegant phenomenological model for cr(r^,xjin 
Eq. (0) which is able to describe all experimental data using only three parameters p7|] . In this model 



o-dipoieir±,x) 



ao I I — exp ^^/^o)^ 



(4) 



with o-Q = 23.03m6, Ql 
description. 



lGey^ xo = 0.0003 and A = 0.288. Figs. | and | show the quahty of this 




Q^ (GeV^) 



FIG. 3. 



Note that even though the above description is impressive, it cannot fix the value of the saturation scale 
from the data which is too constrained by the kinematics ( see Ref. Q for details). 

Thefore, we can conclude that the saturation hypothesis is compatible with all experimental data. However, 
the puzzling situation is that the same data can be described fron a different point of view without a 
saturation scale in the standart DGLAP evolution equation for Q^ > IGeV'^ and the soft phenomenology 
for Q^ < IGeV^. We do not claim that it is a reasonable or smooth parameterization of the data but 
Donachie-Landshoff mixture of soft and hard Pomeron shows that we can produce such a model. 
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FIG. 4. 



Thus in order to fix the saturation scale and to descriniinate between competing models, we need either 
to reach a much smaller values of x ( higher energy) or to use a new target. Realistically, we can conclude 
that 

We need DIS with nuclei to check whether the indications on saturation effect at HERA are really true. 



III. SCALE OF THE SATURATION EFFECT FOR DIS WITH NUCLEI 



A. Asymptotic predictions. 

Let us start with listing the asymptotic predictions of our approach 
Mueller formula for a{r\,x) in Eq. (|) §. 

• At fixed r±_ and at a; ^ 



which is based on the Glauber- 



where Ra is the nucleus radius and h is the surface thickness in the Wood-Saxon nucleon density 



In the same limit 






F^ 



h 
Ra 



OC Q2 ^^. 



The ratio of the diffraction to the total cross sections should depend on energy only weakly ||18| ; 



iif fraction / =^ 



'tot 



il*A) 



Const{W) 



(slowly) -; 



crtot{l*A) 
The energy behaviour of o-''*//''''^*'°"(^*A) is determined by short distances r± « 1/Qs{x; A); 



The high density effects should be stronger in the diffractive channels. 



B. xGa(x,Q^) 

In Fig. 5 we present our calculation of the gluon structure function for different nuclei. Fig.5-d gives a 
glimpse at what we are taking into account in our approach. Figs. 5-a - 5-c show the comparison of our 
calculations, based on the Glauber-Mueller formula, with the solution of the full equation for hdQCD |1J] [^ 
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Fig. 5-a 



Fig.5-b 
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FIG. 5. 

We can derive two conclusions from Fig. 5: (i) the saturation effect is much stronger for a nucler target 
than for a nucleon, and (ii) our model underestimates the value of the effect for Q^ « IGeV^. Unfortunately, 
we have not finished our estimates for F2 for DIS with nuclei. 



'^ Actually, the equation, suggested in Ref. [^ was solved and plotted in Fig. 5 as the asymptotic solution, but this 
equation in the double log approximation coincides with the correct one pM. 



diffractio 

C. R=^i2t 



This ratio shows us how we are close, or how we are far away, from the asymptotic regime since at very 
high energy it should be equal to i. In Fig. 6 we plotted our calculations for this ratio 11^]. One can see that 
the ratio is larger than for the proton target ( see Fig. || ) but it is still smaller than the limiting value of 
i . This is a very encouraging fact for experiment since we do not want to measure a black disc limit which 
is not sensitive to the theoretical approach. In other words, any model or any theoretical approach will give 
the unitarity limit which we call " black disck limit" . A 11 our theoretical QCD prediction are related to the 
form of the transition from pQCD to the "black disc limit" . 
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FIG. 6. 



D. 



F2 slope ( ^2__j_) 



As we have mentioned, gluon satruration leads to a maximum in the Q^ dependence of F2 slope at 
Q^ = Q'^{x;A) fora fixed value of x. Such a maximum has not been seen in the HERA data and has not 
been anticipated in our estimates of the slope. However, the numerical value of the gluon saturation is rather 
for a nucleus target. Figs. and || display the possible experimental effect. In these figures the value of the 
damping factor (D^) for the F2 slope is plotted. D^ is defined in the following way: 



dF,^ 



(x,Q2) 



rflnQ2 



= D''ix,Q^)A 



dF^N:DGLAP^^Q2) 



dlnQ2 



(5) 



where F, 



N-DGLAP 



is the F2 sructure function for a nucleon in the DGLAP approximation. It turns out that 



the value of the effect is sizable for Q < 10 GeV and strongly depends on A. 
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FIG. 8. 



IV. SEARCHING NEW OBSERVABLES. 



F F" 

A. Maxima in ^ and in -^tt 
Ft F^ 



The main idea of this calculation Ig^l is to show that these ratios have maxima at Q^ = Qf,-^^^{x; A) , 
being plotted at fixed x. We want to claim that Q'^^^{x; A) « Qlix; A). Figs. || and |o| show that such a 
suggestion can be right. 
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B. Higher twists in Fl,t and in fPx- 



It is well known that any structure function ( F2{x, Q^) for example ) can be written 



F2ix,Q'') = F^'{x,lnQ') + '-^ F^'^ {x MQ'') + - + 






aP 



F^'ixMQ') 



(6) 



Terms, which are small in terms of Q^ power , called higher twist contributions. The InQ^ dependence of 
the leading and higher twist structure functions ( F2^ and F^^ in Eq. (^ ) is governed by the evolution 
equations. The DGLAP evolution equations 0] give the InQ^ dependence of the leading twist structure 
function ( -Fjf'"'" ) only. Unfortunately, we know only a little about the higher twist contributions. 
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1. We know the evolution equations for all higher twist structure functions [pif ; 

2. We know the behaviour of the higher twist structure functions at low x p2[ . For example, 

3. We know, that higher twist contributions are needed to describe the experimental data p3| |. 

However, it is difficult to estimate the value of the higher twist contributions. Following Ref. |24|, we 
estimate the value of different twist contribution for e A scattering. 

Fig.|ll| shows that the higher twist contributions for nucleus target become smaller than the leading twist 
one only at Q^ > 5GeV'^ even at a; = 10~^. It gives us a hope to treat them theoretically. 

V. CONCLUSIONS. 

1. We have a solid theoretical approach for eA DIS, but we need more experience in numerical solution 
of the non-linear equation specifically for eRHIC kinematic region; 

2. We know pretty well the scale of SC for eA interaction, but we need more systematic study of 
DGLAP evolution for nuclear structure functions and a special investigation whether the initial parton 
distributions could be calculated for nuclear target from the initial parton distributions for proton ; 

3. Our estimates show that we will be able to see the saturation scale in eA DIS at eRHIC being still 
far away from trivial blackening of high energy interaction with nuclei, but we need to check how close 
our model, which we use inn practise, to theoretical estimates; 

4. The ^2 -slope is a very sensitive observable for the saturation scale, but , unfortunaly, we cannot 
expect a qualitatively different behaviou for the saturation models in comparison with others; 

5. Maxima in ratios of Fl/Ft and F£ /F^ give promising tool to extract the value of saturation scale 
Qs{x; A)^ but we need more study on this subject and, in particular, how the initial parton distribution 
for DGLAP evolution could affect our predictions; 

6. eA DIS is very instructive for separation of leading and higher twist contributions, since the fact 
that typical momentum at which these two contributions become of the same order is growing with A. 
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